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A randomly walking quantum particle searches in Grover’s S{y/N) iterations for a 
marked vertex on the complete graph of N vertices by repeatedly querying an oracle 
that flips the amplitude at the marked vertex, scattering by a “coin” flip, and hopping. 
Physically, however, potential energy barriers can hinder the hop and cause the search 
to fail, even when the amplitude of not hopping decreases with N. We correct for these 
errors by interpreting the quantum walk search as an amplitude amplification algorithm 
and modifying the phases applied by the coin flip and oracle such that the amplification 
recovers the ©(v2V) runtime. 
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1 Introduction 

Most quantum algorithms developed to date are based on four general techniques: quantum 
Fourier transforms [e.g., the Deutsch-Jozsa algorithm [T] and Shor’s algorithm [2]), amplitude 
amplification (e.g., Grover’s algorithm and quantum counting i), quantum simulation 
{e.g.^ approximating the Jones polynomial [S] and solving linear equations 0). and quantum 
walks {e.g., element distinctness [7] and NAND evaluation m- Quantum walks are the 
quantum analogues of classical random walks or Markov chains niiniiii], and their crucial 
role in many quantum algorithms has spurred tremendous experimental effort to realize them 
[nmniiiiiiisiiie]. These physical implementations are not ideal, however, and so substantial 
theoretical work has investigated the effects of noise and errors in quantum walk algorithms 

[nmiiiiiiEo]. 

Recently, |20j considered the effects of potential energy barriers hindering a quantum 
particle from searching for a marked vertex on the complete graph of N vertices, as shown in 
Fig. To review, the vertices label computational basis states of an V-dimensional Hilbert 
space, and the d = N — 1 directions from each vertex label a d-dimensional “coin” Hilbert 
space that is necessary to define a non-trivial quantum walk Enna. Then the system \tfj) 
begins in an equal superpositions over both spaces: 

It^o) \^v) G l^c). 
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Fig. 1. (a) The complete graph with N = 6 vertices. A vertex is marked, as indicated by a double 
circle. Identically evolving vertices are identically colored and labeled, (b) Success probability 
as a function of the number of applications of U for search with N = 1024, o = ^ and 

= 0, 0.02z, and 0.04z corresponding to the solid black, dashed red, and dotted green curves, 
respectively. 


where |s„) = \i)/y/N and |sc) = it evolves by repeated applications 

of 

U={aS + /3I){lN<S)Co){Ra<^Id), ( 1 ) 

where a and j3 are the amplitudes of successfully and failing to tunnel through the potential 
barrier, S is the flip-flop shift [53] that causes the particle to hop from one vertex to another 
and then turn around, Cq = 2|sc)(sc| — Id is the Grover diffusion coin [53], and Ra is an 
oracle that flips the sign of the amplitude at the marked vertex. The evolution of the success 
probability as U is repeatedly applied is shown in Fig. [^ Without the potential barrier 
{i.e., 13 = 0), the success probability reaches 1/2 at n^N applications of V [5S]. As /3 
increases, however, the max success probability decreases such that the runtime is 

retained when (3 = 0{1/'/N), and otherwise the algorithm performs no better than classical. 
Thus the hopping errors must not only decrease with N for the search to be fast, but they 
must decrease sufficiently quickly. 

In this paper, we show how to correct for these errors by modifying the phases that the 
coin Co and oracle Ra use, and this recovers the quadratic speedup so long as \f3\ does not 
approach 1. To do this, we first reinterpret the search problem as an amplitude amplification 
algorithm. Then we choose the phases such that the amplification rotates to the marked 
vertex with high probability. Finally, we show that the Q{'/N) runtime is restored. 

2 Amplitude Amplification 

Figure indicates that the system evolves in a 3D subspace spanned by {|a6), \ba), \bb)}, 
where \xy) indicates the equal superposition over the x vertices pointing towards the equal 
superposition over the y vertices [50]. Then in this basis, the operators in Q are 

/O 1 0\ ^ ^ \ /“I 0 o\ 

5=1 0 0 , (/^ 0 Co) = 0 , (i?a 0 /d) = 0 10 

Vo 0 ij ^0 ^ V 0 0 l/ 


( 2 ) 
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Using these, it is straightforward to show that 


IV'-i) 


1 

y/2N-3 




is an eigenvector of Q since it is a 1-eigenvector of S and (trivially) I, and also (—1)- 
eigenvector of the product {1^ (g) Co){Ra ® Id)- Then consider the 2D subspace orthogonal to 
I'i/'-i), which is spanned by 


|s) 




16a) 


N -2 
N-1 


\bb) = 



k) 


1 

71 


(|a6) - |6a)) 



While these are both orthogonal to \tp-i): they are not orthogonal to each other, and we 
express the magnitude of their overlap as the sine of an angle 9: 


sin 9 = 


72(F-1)' 

Now let us find how the operators ([^ act in this subspace. It is straightforward to show that 

S\w) = -\w), S\w-^) = \w^), 


where \w-^) is the state orthogonal to |w), so S' is a reflection through \w). It is also straight¬ 
forward to show that the coin and query together act by 


{In ® Co){Ra ® Id) |s) = |s), 

{lN<E>Co){Ra^Id)\s^)=-\s^), 

so the coin and query together act as a reflection through Thus the search operator U 

Q without errors is the product of two reflections: 

U\a=i,^=o— ^ {In'Si Co){Ra'S i Id) ■ 

Ra ^s-L 

This is simply Grover’s iterate 13 ES], which is illustrated in Fig. with the reflections 
swapped, which does not affect the asymptotic search probability. Thus 7r/40 « ii^fN /2-\/2 
applications of U rotates |s) to |ty) with probability near 1. Since the initial equal superpo¬ 
sition 

\ipo) = (|a6) -h |6a) -h y/N -2\bb)'^ 

is approximately |s) for large N, the system roughly reaches a success probability 1/2 (from 
the |a6) term in \w)) in tt'/N /2-\/2 steps, in agreement with |25j . 
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1^) 



Fig. 2. One application of Grover’s iterate in the 2D subspace spanned by {|s), |n^}} results in a 
rotation of 29. 


Now with potential barriers, the shift/identity term multiplies by a different phase: 

{aS + I31)\w) = —{a — /3)|w) = 

{aS + /3/)|w-‘-) = (a + P)\w-^) = 

where we parameterized a = cos 4>, j3 = i sin </). Since the global phase does not matter, this 
is equivalent to acting by —on \w) and doing nothing to When </) = 0{1/VN), 

this deviation from an exact phase flip is insufficient alter the 0(-\/iV) runtime, but when cj) 
is scales larger, the algorithm is no better than classical |20) . 


3 Correcting Errors 

We now show how to compensate for effective —phase that {aS + /3/) applies to \w) by 
adjusting the coin and query operators, Cg and Ra- In particular, we now use 

C' =(l + e*’')|Se)(Sc|-/d, 


for some phase t] that we must determine. Then the coin flip and oracle together are 


[In 0 C'g)iK ^ Id) 



0 

N-1 _ 

(l+e”')VJV^ 

N-1 


0 

(l+e‘'')V7V=2 

N-1 

{N-2)e'^-l 

N-1 


Since |V'-i) is still a (—l)-eigenvector of this, we again consider the 2D subspace spanned by 
{|s), |w)}. It is straightforward to show that 


=-|sy. 

But since the global phase does not matter, this is equivalent to multiplying the phase of |s-*-) 
by —6“*’* and doing nothing to |s). 

Thus the action of the modified search operator is 


U'= {aS + f3I) 


QJJ 


— on |s-L) 
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Fig. 3. Success probability as a function of the number of applications of U' (search corrected 
for potential barriers) for search with N = 1024, ol = — |/3p, and /3 = 0, 0.42, and 0.82 

corresponding to the solid black, dashed red, and dotted green curves, respectively. 


H 0 yer derived a condition on these two phases if and only if amplitude amplification rotates 
the state \s) to \w) with certainty |27) . which in our notation is 

tan ^ 2 ^ '^ ~ (2) d ~ 2sin^ 9 ). 


Solving for 77 , which tells us what phase the adjusted coin Cq and query R'^ should use, we 
get 

77 = —2 tan“^ 

With this choice of 77 , Fig. shows the success probability as we repeatedly apply U' with 
N = 1024 and /3 = 0, OAi, and O.Si. Even with these strong potential barriers, our correction 
still causes the system to rotate from {-ipo) « |s) to jic) with probability near 1 , as expected, 
which results in a success probability of 1/2 from the \ab) term in | 7 c). 

As (3 increases, however, the algorithm slows down. We can find the runtime by deter¬ 
mining the angle of rotation a (which was 29 for the barrier-free case in Fig. using Hpyer’s 
result m- sincT = \ {w\U'\tjjo)\: where \tpo) = |si,)( 8 )|sc) is the initial equal superposition state. 
Evaluating this inner product with rj = —2(j), 


tan (</) 


N-1 

N-2 


-2(j). 


{w\U'\tpo) 


(1 g-2*0)g-z0 

y/m 


The magnitude of this is 


sinCT = \{w\U'\tpo)\ = 


1 -I- cos 2(/) 

N ■ 


We want to rotate by a total amount roughly tt/2, so if the number of applications of U' is 
t*, then we want ut* = 7r/2, or 


TT 


2 sin 


-1 ( / l+cos 20 ^ 2 -v^l -|- cos 2(j) 

N 


y/N 
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for large N. This runtime agrees with Fig. e.g., with N = 1024 and (3 = 0.8i ^ (j) = 
sin“^(0.8), we get = 59, which corresponds to the peak in success probability. 

Using this runtime formula, we now show that the quadratic quantum speedup is recovered 
when (f) does not approach n/2. First, when (p scales less than a constant, then cos 2(/) « 1 for 
large N, so the runtime is simply the barrier-free tt'/NI 2\f2. Furthermore, when (p scales as 
a constant (less than 7r/2 so the particle still hops), then while the runtime is slower, it still 
has the same Q{'/N) scaling. Thus we have corrected for potential barriers by changing the 
phases of the coin flip Cq and query Ra in such a way that amplitude amplification recovers 
the full quadratic speedup, so long as the potential barriers get no worse as the problem grows 
(and assuming ^ ^ 7r/2 so that the particle still hops). 

In the other extreme when p approaches 7 r/ 2 , we have p = 12 — 6 with i5 —>■ 0. In this 

case, the particle mostly stays put, and one can consider it an “error” for the particle to 
tunnel to another vertex. Then our formula for yields a runtime of tt'/N /2-\/2(5. So when <5 
scales less than the algorithm still has some improvement over classical. Thus rather 

than quantifying how small potential barriers must be to allow fast search, this quantifies how 
high potential barriers must be to stop the search. Finally, when (5 = 0 so that p = 7 r/ 2 , the 
particle does not hop at all, and is infinite, as expected, since the state is never rotated to 
\w). 

4 Continuous-Time Quantum Walk 

For completeness, we briefly state the correction for continuous-time quantum walks with 
potential barriers, which was given in |20j . If the potential barriers lower the transition 
amplitude of the particle by e, then the jumping rate 7 can be adjusted from its barrier-free 
value of \/N to 

Af(l-e)’ 

and the walk searches with probability I in time /2 as if there were no potential barriers. 

5 Conclusion 

In the standard formulation of quantum walk search, even small potential barriers affecting 
the particle’s hop can eliminate the quadratic quantum speedup. We have shown, however, 
that these errors can be corrected for by modifying the phases of the coin and query operators. 
By interpreting the quantum walk as an amplitude amplification algorithm, we find the precise 
phases to use and show that the quadratic speedup is restored. 

This approach to correcting errors from potential barriers should be applicable to a variety 
of quantum walk algorithms, so long as the evolution approximately occurs in a 2D subspace 
admitting rotations as in Grover’s algorithm in Fig. such as element distinctness [7]. For 
algorithms where this does not hold, how to correct for potential barriers remains an open 
question, as does the effect of the potential barriers themselves. 

Throughout this paper, we have assumed that the potential barriers are fixed and identical 
between every pair of vertices. Physically, however, they may be non-uniform or fluctuate 
randomly. The effects of these generalizations are open questions, as well as how to correct 
for them. They would also likely break the symmetry that makes this work tractable. 
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